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Abstract 



This paper is a short introduction to orthogonal polynomials, both the general theory 
and some special classes. It ends with some remarks about the usage of computer algebra for 
this theory. The paper will appear as a chapter in the book "Computer Algebra in Quantum 
Field Theory: Integration, Summation and Special Functions" , Springer- Ver lag. 

1 Introduction 

This paper is a short introduction to orthogonal polynomials, both the general theory and 
some special classes. After the definition and first examples in Section 2, important but mainly 
elementary aspects of the general theory associated with the three-term recurrence relation 
are treated in Section 3. Sections 4, 6 and 7 discuss special classes of orthogonal polynomials, 
interrupted by Section 5 about Gauss quadrature. Section 8 collects some more advanced results 
in the general theory of orthogonal polynomials. Finally Section 9 discusses the role of computer 
algebra in the theory of (special) orthogonal polynomials. 

Everything treated here is well-known from the literature. I mention a few books which 
can be recommended for more detailed study. A great classical introduction to orthogonal 
polynomials, both the general theory and the special polynomials, is Szegd [Mj- A very readable 
textbook, in particular for the general theory, is Chihara [5]- As a textbook emphasizing the 
special theory I recommend Andrews, Askey & Roy [2]. Very good is also Ismail [10], but more 
focusing on the g-case. Two recent compendia of formulas for special orthogonal polynomials 
are Giver et al. [ISt Gh. 18] and Koekoek, Lesky & Swarttouw pi Gh. 9 and 14]. 

2 Definition of orthogonal polynomials and first examples 

Let V be the real vector space of all polynomials in one variable with real coefficients. Assume on 
V a (positive definite) inner product (/,(?) (7,5 G V). Orthogonalize the sequence pf monomials 
1, X, , . . . with respect to the inner product (Gram-Schmidt). This results into the sequence 
Pq{x)^Pi{x),P2{x), ... of polynomials in x. So po{x) = 1 and, if pq{x)^pi{x)^ . . . ,p„_i(x) are 
already produced and mutually orthogonal, then 

Pn{x):=x - > 7 rPfc(a;). 
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Indeed, p„(x) is a linear combination of 1 and 



n-l 



{Pn,Pj) = {x'^,Pj) - ^ {Pk,Pj) = {x'^,Pj) - (Pj'Pj) = (i = 0, 1, . . . ,n - 1). 

k=0 ^Pk'Pk) \Pj^Pj/ 

Throughout we will use the constants hn and associated with the orthogonal system: 

{PmPn) = hn, ^^(x) = fc^x" + polynomial of lower degree. (2-1) 

The Pn are unique up to a nonzero constant real factor. We may take them, for instance, 
orthonormal {hn = 1; this determines p„ uniquely if also kn > 0) or monic {kn = 1). 
In general we want 

{xf,g) = {f,xg). 

This is true, for instance, if 

pb °° 
{fi9)-= f{x)g{x)w{x)dx or {f,g)-='^f{xj)g{xj)wj 

i=o 

for a weight function w{x) > or for weights Wj > 0, respectively. These are special cases of an 
inner product 

{f,g) ■■= I f{x)g{x)dfi{x), (2.2) 



where ^ is a (positive) measure M, namely the cases d^{x) = w{x) dx on an interval /, and 
= H'jLi , respectively. 

A measure /U on K can also be thought as a non- decreasing function fx on M. Then 

/ f{x)dfi{x)= / f{x) djl{x) = lim / f{x)djl{x) 

can be considered as a Riemann-Stieltjes integral. The measure fi has in x a mass point of mass 
c > if the non-decreasing function Jl has a jump c at x, i.e., if lim^j^g {jj-{x+5)—'jl{x—6)) = c > 0. 
The number of mass points is countable. More generally, the support of the measure fi consists 
of all X e M such that pt(x + 5) — Jl{x — (5) > for all 5 > 0. This set supp(^) is always a closed 
subset of M. 

In the most general case let /i be a (positive) measure on M (of infinite support, i.e., not 
M = -LU i'. ) ™* / "M.) < oo for all = 0, 1, 2 A system . . .) 

JR 

obtained by orthogonalization of {1, x, x^, . . .} with respect to the inner product ()2.2p is called 
a system of orthogonal polynomials with respect to the orthogonality measure ii. 
Here follow some first examples of explicit orthogonal polynomials. 

• Legendre polynomials Pn{x), orthogonal on [—1,1] with respect to the weight function 1. 
Normalized by Pni^) = 1- 
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• Hermite polynomials Hn(x), orthogonal on (—00,00) with respect to the weight function 

2 

. Normahzed by /c„ = 2". 

• Charlier polynomials c„(x, a), orthogonal on the points x = 0, 1, 2, . . . with respect to the 
weights a^/x! (a > 0). Normalized by c„(0;a) = 1. 

The hn (see (j2.ip ) can be computed for these examples: 

Pmix) Pn{x) dx = — -(5m,n, 7r"2 / i?^ (x) ff„ (x) e"^ = 2"n! (5m,,n , 

-1 zn + i 

00 ™ 

e"" ^ Cmix, a) c„(x, a) — = a~"n! 5m,n ■ 

3 Three-term recurrence relation and some consequences 
3.1 Three-term recurrence relation 

The following theorem is fundamental for the general theory of orthogonal polynomials. 
Theorem 3.1. Orthogonal polynomials Pn satisfy 

xpn{x) = anPn+l{x) + bnPn{x) + C„p„_i(x) (n > 0), 

xPq{x) = aopi(x) + bQPQ{x) 
with an,bn,Cn real constants and a„c„+i > 0. Also an = , — ■ 

kn+l hn+1 hn 

Moreover (Favard theorem), if polynomials pn of degree n (n = 0, 1, 2, . . .) satisfy (13. ip with 
an,bn,Cn real constants and a„c„+i > then there exists a {positive) measure on R such that 
the polynomials pn are orthogonal with respect to fi. 

The proof of the first part is easy. Indeed, xpn{x) = X^fcio ^kPk{x), and if A; < n — 2 then 
{xpn,Pk) = {Pn,xpk) = 0, hencc Uk = 0. Furthermore, 

^ _ {xpn+l,Pn) _ {xPn,Pn+l) _ {xPn,Pn+l) hn+1 _ ^ hn+1 
"'''^ {Pn,Pn) hn {Pn+l,Pn+l) hn " hn 

Hence a„c„+i = a^ hn+i/hn > 0. Hence c„+i//i„+i = a„//i„. 

The proof of the second part is much deeper (see Cihara [5] Ch. 2]). 
Remarks 



(3.2) 



1. For orthonormal polynomials the recurrence relation (j3.ip becomes 

XPn{x) = anPn+l{x) + bnPn{x) + a„_ip„_i(x) (n > 0), 

xpo(x) = aopi(x) + bQPQ{x), 
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and for monic orthogonal polynomials 

XPn{x) = Pn+l{x) + bnPn{x) + CnPn-lix) (n > 0), 

xpo{x) = pi{x) + bopo{x), 

with Cji = hn/hn-i > in (I3.3p . If orthonormal polynomials Pn satisfy (I3.2p then the 
corresponding monic polynomials k~^pn satisfy (|3.3|) with c„ = a^_i • 

2. If the orthogonality measure is even {fi{E) = fi(—E)) then Pn{—x) = {—l)"'Pn{x), hence 
bn = 0, so xpn{x) = a„j)„+i(x) + CnPn-iix)- Examples of orthogonal polynomials with 
even orthogonality measure are the Legendre and Hermite polynomials. 

3. The recurrence relation ()3.ip determines the polynomials pn uniquely (up to a constant 
factor because of the choice of the constant po). 

4. The orthogonality measure fi for a system of orthogonal polynomials may not be unique 
(up to a constant positive factor). See Example 13.21 

5. If /i is unique then V is dense in L^(/u). See Shohat &: Tamarkin [2H Theorem 2.14]. 

6. If there is an orthogonality measure fi with bounded support then fi is unique. See Chihara 
[3 Ch. 2, Theorem 5.6]. 



3.2 Moments 

The moment functional M : p ^ (p, 1) : "P — )• M associated with an orthogonality measure is 
already determined by the rule M{pn) = {pn, 1) = for n > 0. Hence M is determined (up to 
a constant factor) by the system of orthogonal polynomials pn, independent of the choice of the 
orthogonality measure, and hence M is also determined by p. ID . The same is true for the inner 
product {f,g) = {fg,l) on V. 

The moment functional M is also determined by the moments fin ■= {x"", 1) (n = 0,l,2,...). 
The condition o„Cn+i > is equivalent to positive definiteness of the moments, stated as 

An := det{fii+j)lj^Q > (n = 0, 1, 2, . . .). 

For given moments fin and corresponding orthogonal polynomials pn a positive measure fj, is an 
orthogonality measure for the pn iff is a solution of the (Hamburger) moment problem 

[ x^dfi{x)=fin (n = 0,1,2,...). (3.4) 

Uniqueness of the orthogonality measure is equivalent to uniqueness of the moment problem. 

Example 3.2 (non-unique orthogonality measure). The following goes back to Stieltjes \1?>\ 
§56]. In the easily verified formula 

/oo 
e-"'(l + C7sin(27rn)) (in = vr^/^ 
-oo 
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make a transformation of integration variable u = logx — ^(n + 1) and take — 1 < C < 1. Then 

POO 

/ x"(l + Csin(27rlogx))e-i°s'^dx = e3"("+2)_ (35) 

^0 



TT 2 e 4 



Thus a one-parameter family of measures yields moments which are independent of C. The 
corresponding orthogonal polynomials pn are a special case of the Stieltjes- Wigert polynomials 
[I2I §14.27]: Pn{x) = Sniq'^x^q) with q = e~2^ see Christiansen |6l p. 223]. 
It is also elementary to show that 



00 



\k=—OD / k=—oo 

which means that the same moments, up to a constant factor, as in (13. Sh are obtained with the 

^exp(-ifc)- 



measure Er=-oo e ^^'^+'^''5exp(-i/ 



3.3 Christoffel-Darboux formula 

Let Vn be the space of polynomials of degree < n. Let {pn} be a system of orthogonal polyno- 
mials with respect to the measure /i. The Christoffel-Darboux kernel is defined by 

K.ix,y):=±P-i^^. (3.7) 
j=0 ^ 

Then 

(n„/)(x) := / Kn{x,y)f{y)dfi{y) 

defines an orthogonal projection n„: V ^ Vn ■ Indeed, if /(y) = J2T=o '^kPkiu) (finite sum) 
then 

n 00 „ n 

(n„/)(x) = J^Pi'^^^^ir / Pj(y)Pkiy)dKy) = Y(^jPjix)- 

The Christoffel-Darboux formula for Kn{x,y) given by ()3.7p is as follows. 



Vj{x)Pj{y) 



hj 



( kn Pn+l{x)pn{y) - Pn{x)pn+l{y) , ^ . ,^ 

hnkn+i x-y 
k 

{Pn+lix)Pn{x) - p'n{x)pn+l{x)) {x = y). (3.9) 



For the proof of (|3.8p note that 



xpj{x) = ajpj+i{x) + bjpj{x) + CjPj-i{x), 

ypj{y) = ajPj+i{y) + bjPjiy) + CjPj-i{y). 
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Hence 



hj hj hj 



Use that Cj/hj = aj^i/hj^i. Sum from j = to n. Use that a„ = kn/kn+i. This yields (j3.8p . 
For dSl]) let y ^ in 1^^. 



3.4 Zeros of orthogonal polynomials 

Theorem 3.3. Let pn he an orthogonal polynomial of degree n. Let /x have support within the 
closure of the interval {a,b). Then pn has n distinct zeros on {a,b). Furthermore, the zeros of 
Pn and Pn+i alternate. 

Proof For the proof of the first part suppose pn has precisely k < n sign changes on (a, b) at 
xi,X2, ■ ■ ■ ,Xk- Hence, after possibly multiplying p„ by —1, we have pn{x){x — xi) . . . {x — Xk) > 
on [a, b]. Hence Jl^pn{x){x — xi) . . . (x — x^) dfi{x) > 0. But by orthogonality we have 
JaPn{x){x — xi) . . . (x — Xfc) dfi{x) = 0. Contradiction. □ 
For the proof of the second part use (j3.9p : If kn, kn+i > then 

Pn+lix)Pnix) - p'^{x)pn+lix) = " ^ ^7 > 0. 

j=0 

Hence, if y, z are two successive zeros of Pn+i then 

P'n+l^y)Pn{y) > 0, p'n+l{z)pn{z) > 0. 

Since p'^_^_l{y) and have opposite signs, Pn{y) andp„(z) must have opposite signs. Hence 

Pn must have a zero in the interval {y, z). □ 

3.5 Kernel polynomials 

Recall the Christoffel-Darboux formula ()3.8p . Suppose the orthogonality measure p, has support 
within (—00,6] and fix y > 6. Then 

Kn{x, y) x^ {y - x) dp{x) = y^{y - y) = {k <n). 

Hence x i— )• qn{x) := Kn{x,y) is an orthogonal polynomial of degree n on (—00,6] with respect 
to the measure (y — x) dp{x). Hence 



/ \ I \ Pn{y) I \ 

qn{x) - qn-l{x) = — Pn{x), 

hn 



hnkn+1 



Pn{y)Pn+l{x) - Pn+l{y)Pn{x) = 7 {x - y)qn{x). 

Kn 

The orthogonal polynomials qn are called kernel polynomials. Of course, they depend on the 
choice of p and of y. 
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4 Very classical orthogonal polynomials 



These are the Jacobi, Laguerre and Hermite polynomials. They are usually called classical 
orthogonal polynomials, but I prefer to call them very classical and to consider all polynomials 
in the (g-)Askey scheme (see Sections 6 and 7) as classical. 
We will need hypergeometric series [2, Ch. 2]: 

fai,...,ar \ _ ^ {ai)k...{ar)k f... 

'•^\bs,...,bs'V'~t,{h)k...{bs)k kr 

where (a)^ := a{a + 1) ... (a + A; — 1) for k = 1,2, . . . and (a)o := 1 is the shifted factorial. If 
one of the upper parameters in ()4.ip equals a non-positive integer — n then the series terminates 
after the term with k = n. 



4.1 Jacobi polynomials 

Jacobi polynomials Pn"'^^ [iSl §9.8] are orthogonal on (—1,1) with respect to the weight function 
w{x) := (1 — x)"(l + (cK,/3 > —1) and they are normalized by P^"'^^(l) = (a + l)„/n! . They 
can be expressed as terminating Gauss hypergeometric series: 

P,.,«(,) = (2i±iV " ; ;| " + S i(i - .)) 

^ A (n + g + /3 + l)fc {a + k + l)„_fc {x - 1)^ 

^ k\in-k)\ 2k ■ ^ ■ ' 

They satisfy (because of the orthogonality property) the symmetry P^'^^ {—x) = (—1)" Pi^'°'\x). 
Thus we conclude (much easier than by manipulation of the hypergeometric series) that 

/_n,n + a + /3 + l \ (-l)"(/3 + l), /-n,n + a + ^ + l 
'"^'K a + 1 '"j- (a + 1). '"^"K /3 + 1 

For Pn{x) := P^'^\x) there is the second order differential equation 

(1 - x^)p'^{x) + {13 -a- {a + 13 + 2)x)p'^{x) = -n{n + q + /3 + l)p„(x). (4.3) 
This can be split up by the shift operator relations 

^ Pi'-^'X^) = \{n + a + P + l)PtV''^'\x), (4.4) 
(1 - ^ Ptf''^'\x) + {(3-a-{a + P + 2)x)PtV''^'\x) 

= (1 - x)-"(l + x)-'^ ^ ((1 - xr+\l + xf+^PtX^'^+^\x)) = -2nPt'^\x). (4.5) 
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Note that the operator d/dx acting at the left-hand side of (|4.4p raises the parameters and 
lowers the degree of the Jacobi polynomial, while the operator acting at the left-hand side of 
(|4.5p lowers the parameters and raises the degree. Iteration of (|4.5|) gives the Rodrigues formula 



Special cases 



Gegenbauer or ultraspherical polynomials {a = (3 = X—^): C^(x) := — — ^'"^ ^^(^)- 



• Legendre polynomials (a = /3 = 0): Pn{x) := Pji''^\x). 



• Chebyshev polynomials (a = /3 = it^): 

^ f Q\ f a\ ^' T3(-h~l)t Q\ TT I a\ sin(n + l)6' (2)„ (1,1) 

Tn[cos9) :=cos{n9) = —r—Pn (cosfc*), C/„(cos6') := — = Pn [cosU). 

(^)n sm6' (|)„ 

Quadratic transformations Since P2n'°'\^) even polynomial of degree 2n in x, it is 

also a polynomial p„(2a;^ — 1) of degree n in 2x^ — 1. For m / n we have 

= ^ p„(2y2 _ l)p„(2y2 - 1) (1 - y^)" = const, j ^p,n{x)pn{x) (1 - + x)"^ dx. 

Hence 



Similarly, 





_pi"'-^\2x2- 


1) 




Pi"'~~'\l) 




-f^in+l (^) 


_ xPt'^\2x^ - 


1) 









Theorem 4.1. [5| Ch. 1, §8] Let {pn} be a system of orthogonal polynomials with respect to an 
even weight function w on M. Then there are systems {qn} and {r„,} of orthogonal polynomials 
on [0,oo) with respect to weight functions x ^ x^2u){x'2) and x i— )■ x'2w{x2), respectively, such 
thatp2n{x) = Qnix'^) and p2n+i{x) = xr„(x2). 

4.2 Electrostatic interpretation of zeros 

Let Pn{x) := const. Pn^ ^'^"^ '^\x) = {x — xi){x — X2) . . .{x — x^) be monic Jacobi polynomials 
(p,?>0). By (gJl) 

(1 - x'^)p'l^{x) + 2{q-p - {p + q)x)p'n{x) = -n{n + 2p + 2q- l)p„(x). 
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Hence 



= 1- Xk)pl{xk) +2{q-p-{p + q)xk)p'n{xk) 



1 Pn{Xk) , P , ^ - = ^ 1 , P 



2 ' - 1 ' Xk + 1 /-^^ Xk-Xj ' Xfc - 1 ' Xk + l' 

This can be reformulated as 



(Vy)(xi,...,rc„) =0, 

where 

i<j j j 

is the logarithmic potential obtained from n + 2 charges q,l, . . . ,l,p at successive points —1 < 
j/i < . . . < y„ < 1. It achieves a minimum at the zeros goes back to 



Stieltjes [22]. 

4.3 Laguerre polynomials 

Laguerre polynomials L" |12t §9.12] are orthogonal on [0,oo) with respect to the weight function 
w{x) := x"e~^ (a > —1). They are normalized by -Z^^(O) = (a + l)n/'n\ . They can be expressed 
in terms of terminating confluent hypergeometric functions by 

Lnix) = ^-^^ iFi ( ~\;x)=y: ^""uu ^ nT' (4-6) 
^ n! \a + l J kl{n-k)\ ^ ' ^ ^ 

For Pn{x) := there is the second order differential equation 

xPn{x) + (a + 1 - x)p'j^{x) = -npn{x). 
This can be split up by the shift operator relations 

— = — L„^i(x), 

and 

X ^ Llt\{x) + (a + 1 - x)Llt\{x) = x-e^ [x'^^^ Llt\{x)) =nLl{x). (4.7) 
Iteration of (j4.7|) gives the Rodrigues formula 

-a x f A 
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4.4 Hermite polynomials 

Hermite polynomials Hn |12^ §9.15] are orthogonal with respect to the weight function w{x) := 

2 

on M and they are normalized by Hn = 2"rE" + • • • . They have the explicit expression 

gn(x) = n!^ ^ . (4.8) 

^ j!(n-2j)! 

There is the second order differential equation 

H'^{x) - 2xH'n{x) = -2nHn{x). 
This can be split up by the shift operator relations 

H'n{x) = 2nHn-i{x), K_^{x) - 2xHn-i{x) = e"' ^ (e-"'//„_i(x)) = -F„(x). (4.9) 

Iteration of the last equality in (j4.9p gives the Rodrigues formula 

/ J \ " , 



Hn{x) = {-lTe^ (-) [e- 
4.5 General method to derive the standard formulas 

The previous formulas can be derived by the following general method. Let (a, h) be an open 
interval and let w^wi > be strictly positive C^-functions on (a, 6). Let {pn} and be 
systems of monic orthogonal polynomials on (a, h) with respect to the weight function w resp. 
wi. Then under suitable boundary assumptions for w and wi we have 

p'n{x) qm-i{x) wi{x) dx = - J pn{x)w{xy^ —{wi{x) qm-i{x)) w{x)dx. 

Suppose that for certain a„ / : 

w{x)~^ — {wi{x)x'^~^) = —anx'^ + polynomial of degree < n. 
dx ' 

Then we easily derive a pair of first order differentiation formulas connecting {pn} and {(/n}; 
an eigenvalue equation for p„ involving a second order differential operator, and a formula 
connecting the quadratic norms for p„ and qn-\ ■ 

Pn{x) = nqn-i{x), w{xy^ {wi{x) qn-i{x)) = -anPn{x), 

dx 

■w{xy^ (wiix) Pn{x)) = -nanPnix), 

/•b pb 

n I qn-iix)"^ wi{x) dx = Un / Pn{x)'^ w{x) dx. 
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In particular, if we work with monic Jacobi polynomials Pri'^\ then (a, b) = ( — 1, 1), w{x) = 
{l-xr{l + xf,Pn{x)=pt^\x), wi{x) = {l-xr+^{l+xf+\ qm{x)=pt^^'^^^\x). Then 
ttri = {n + a + (5 + 1). Hence 

^pt'\x) = n^:X"''^'\x), (4.10) 
l-x')^ + {(3-a-{a + l3 + 2)x))pi"+/'''+'^(x) = -{n + a + l3 + l)pt^\x). (4.11) 
For x = l M.llj) yields 

Iteration gives 

^" (n + a + /3 + l)„- ^^-''^ 

So for p„ = const. pn"''''' = k^x^ + • • • we know pn{l)/kn, independent of the normalization. 

The hypergeometric series representation of Jacobi polynomials is next obtained from (|4.10p 
by Taylor expansion: 

P^^^(.) _A(x-l)^ /d^ 



_^ {x-lf nl Prr""'^'a) _ p/-n,n + a + /3 + l , \ 
k\ {n-k)\ ^(-./^)m ' \ a + 1 



x=l 

(o+fc,/3+fc) 
n~fc 

The quadratic norm /i„ can be obtained by iteration of 
1 n 



X Ti + a + p + l y 

So for pn = const, pn"'''^ = + • • • we know hn/k'^, independent of the normalization. 



4.6 Characterization theorems 

Up to a constant factors and up to transformations x — )■ ax + 6 of the argument the very classical 
orthogonal polynomials (Jacobi, Laguerre and Hermite) are uniquely determined as orthogonal 
polynomials Pn satisfying any of the following three criteria. (In fact there are more ways to 
characterize these polynomials, see Al-Salam JT\.) 

• {Bochner theorem) The pn are eigenfunctions of a second order differential operator. 

• The polynomials p'^+i are again orthogonal polynomials. 
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• The polynomials are orthogonal with respect to a positive weight function w on an 
open interval / and there is a polynomial X such that the Rodrigues formula holds on /: 



Pn{x) = const. w(x) 




4.7 Limit results 

The very classical orthogonal polynomials are connected to each other by limit relations. We 
give these limits below for the monic versions pC'^\i'^,hn of these polynomials, and on each 
line we give also the corresponding limit for the weight functions: 

^n/2^(a,")(^/c,i/2) ^ /,^(^)^ (1 " X Va)° ^ e"^'' , « ^ oo, (4.13) 

{-I3/2Yp^^^^\1-2x/P) ^ C{x), x"(l-x//3)'^ ^x°e-^\ /3 ^ oo, (4.14) 

(2a)~"/2^^((2a)^/2^ + a) ^ /i„(x), {I + {2/af/^x)^e-'^'^°'^''''' ^ e-''\ a^oo. (4.15) 

The limits of the orthogonal polynomials in ()4.13p and (I4.14|] immediately follow from (14. 2p . 
(|4.6p and (|4.8p . For various ways to prove ()4.15p see [13, section 2]. 



5 Gauss quadrature 

Let be given n real points xi < X2 < ■■■ < x„,. Put Pn{x) := {x — xi) . . . {x — x„). For 
fc = 1, . . . , n let Zfc be the unique polynomial of degree < n such that lk{xj) = 6k j {j = 1, . . . , n). 
This polynomial, called the Lagrange interpolation polynomial, equals 



Uj-j^Xk- Xj) {X-Xk)p'n{xk) 

For all polynomials r of degree < n we have 

n 

r{x) = ^r{xk)lk{x). 

k=l 

Theorem 5.1 (Gauss quadrature). Let pn be an orthogonal polynomial with respect to /i and 
let the Ik be the Lagrange interpolation polynomials associated with the zeros xi, . . . ,x„ of Pn-- 
Put 

Afc := / lk{x)dii{x). 



Then 

\2 



lk{x) dfi{x) > 
and for all polynomials of degree < 2n — 1 we have 

„ n 

/ f{x)dfi{x) = y^Xkf{xk). (5.1) 



k=l 
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Proof Let / be a polynomial of degree < 2n — 1. Then for certain polynomials q and r of 
degree < n — 1 we have f{x) = q{x)pn{x) + r{x). Hence f{xk) = r{xk) and 



Also 



„ n „ n n 

)= / r{x)d^{x) =^r{xk) / lk{x) dii{x) = ^\kr{xk) = ^\kf{xk)- 
■'^ k=i ''^ k=i k=i 

n „ 

^k = y] Xj lk{xjf = / lk{xf dfi{x) > 0. □ 
,=1 



3 

Prom (j5.ip we see in particular that, for i, j < n — 1, 



„ n 

= / Piix)Pj{x)dfi{x) =y2^kPiixk)Pjixk)■ 
■^^ k=l 

Thus the finite system po,Pi, ■ ■ ■ ,Pn-i forms a set of orthogonal polynomials on the finite set 
{xi, . . . , Xn} of the n zeros of pn with respect to the weights \k and with quadaratic norms hj. 
All information about this system is already contained in the finite system of recurrence relations 

xpj{x) = ajPj+i{x) + bjpj{x) + CjPj-i{x) (j = 0, 1, . . . ,n - 1) 

with ttjCj+i > (j = 0, 1, . . . , n — 2). In particular, the Xk are obtained up to a constant factor 
by solving the system 

n 

^^XkPjixk) = (j = 1, . . . ,n - 1). 



k=l 



6 Askey scheme 

As an example of a finite system of orthogonal polynomials as described at the end of the previous 
section, consider orthogonal polynomials po,pi, . . . ,pN on the zeros 0, 1, . . . , of the polynomial 
PN-<ri{x) ■= x{x — 1) ... (x — N) with respect to nice explicit weights Wx (x = 0, 1, . . . , A) like: 

• Wx '■= — p)^~^ (0 < p < 1). Then the Pn are the Krawtchouk polynomials 

f An 77^ f-n,-x 1\ ^ {-n)k{-x)k 1 

• Wx ■= ^'^ -^r~' — ^~Tr~ (cK)/^ > Then the pn are the Hahn polynomials 

x! (A*— xj! 

^ / o AT\ 77 f-n,n + a + /3 + l,-x 
\ a + 1, —A 
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Hahn polynomials are discrete versions of Jacobi polynomials: 

r, (AT R AT\ „ f -n,n + a + P + 1,-Nx 
Qn[Nx;a,(3,N) = 3F2( , . 

V a + 1, —N 



2F1 ^ ^ ; X = const. Pr^\l - 2x) 

\ a + 1 / 



and 



Qm{Nx;a,P,N)Qn{Nx;a,P,N)wNx 

xe{oi,^,...,i} 



const, r - 2x)Pi"'^)(l - 2a;) x°(l 

Jo 



x)^ dx. 



Jacobi and Krawtchouk polynomials are different ways of looking at the matrix elements 
of the irreducible representations of SU(2), see [16]. The 3j coefficients or Clebsch-Gordan 
coefficients for SU(2) can be expressed as Hahn polynomials, see for instance |15j. 

While we saw that the Jacobi, Laguerre and Hermite polynomials are eigenfunctions of a 
second order differential operator, 

A{x)p'^{x) + B{x)p'^{x) + C{x)pn{x) = XnPnix), (6.1) 

the Hahn and Krawtchouk polynomials are examples of orthogonal polynomials pn on {0, 1, . . . , N} 
which are eigenfunctions of a second order difference operator, 

A{x)pnix - 1) + B{x)pnix) + C(x)p„(x + 1) = A„p„(x). (6.2) 

If we also allow orthogonal polynomials on the infinite set {0, 1,2,...} then Meixner polyno- 
mials Mn{x; 13, c) and Charlier polynomials Cn{x]a) appear. Here 

M„(x;/3,c):= 2fJ~\~''-1- nj,:=^Mc^^ 
VP c J x! 

Cn{x; a) := 2Fo{-n, -x; ; -a~^), := a^/x! . 

If we also include orthogonal polynomials which are eigenfunctions of a second order operator 
as follows, 

A{x)pn{x + i) + B{x)pn{x) + C{x)pn{x - i) = XnPn{x), (6.3) 

then we have collected all families of orthogonal polynomials which belong to the Hahn class. 
Similarly, with an eigenvalue equation of the form 

A{x)pniq{x + 1)) + B{x)pniq{x)) + C{x)pn{q{x - 1)) = \nPn{q{x)), (6.4) 

where g is a fixed polynomial of second degree, we obtain the orthogonal polynomials on a 
quadratic lattice. All orthogonal polynomials satisfying an equation of the form (j6.ip - (l6.4p 
have been classified. There are only 13 families, depending on at most four parameters, and 
all expressible as hypergeometric functions, 4F3 in the most complicated case. They can be 
arranged hierarchically according to limit transitions denoted by arrows. This is the famous 
Askey scheme, see for instance [17, Fig.l]. 
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7 The g-case 



On top of the Askey-scheme is lying the g-Askey scheme |12^ beginning of Ch. 14] , from which 
there are also arrows to the Askey scheme as g — t- 1. We take always < g < 1 and let g t 1 
for the limit to the classical case. Some typical examples of g-analogues of classical concepts are 
(see Gasper &: Rahman [3]): 

1 - 

• g- number: [a]„ := > a 

I — q 

fa. \ 

• g-shifted factorial: (a;q)n := TT(1 — aq'') (also for n = oo), - — ^ — > (a)k- 



g-hypergeometric series: 



ai,...,as+i \ ^ (ai;g)fc . . . (as+i;g)fc 



bi,...,bs ' ' J ' ^ {bi;q)k...{bs]q)k iq;q)k' 
g"i,...,g"''+i \ ^ / ai,...,as+i 



'^''^'^ q'^,...,q'^ ;g,zj^.-MF.^ ^^^^^^^^^ 



• g-derivative: {Dgf){x) := — r > f'{x). 

(1 - q)x 

. g-integral: / /(x) := (1 - g) V /(g'^) g'^ ^ / 

^0 ._n ^0 



/(x) dx. 



The g-case allows more symmetry which may be broken when taking limits for g to 1. In 
the elliptic case [9, Ch. 11] lying above the g-case there is even more symmetry. 

On the highest level in the g- Askey scheme are the Askey- Wilson polynomials [3]. They are 
given by 

{ab;q)n{ac;q)n{ad;q)n , f q~'^,q''~'^abcd,ae''^ ,ae~'^ \ 

Pn{cos 6; a, b,c,d\ q) := 4</>3 , , ■,q,q] , 

a"- \ ao, ac, ad J 

and they are symmetric in the parameters a, 6, c, d. For suitable restrictions on the parameters 
they are orthogonal with respect to an explicit weight function on (—1, 1). In the special case 
a = —c = b = —d = {qf3)2 we get the continuous q-ultraspherical polynomials |12l §14.10.1]. 
They satisfy the orthogonality relation 



Cm (cos 6*; /3 I g)C„(cos6';/3 | g) 



(e2^^;g)oc ' 



(/3e2'^;g)c 



de = (m^n) 
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and they have the generating function 

^ 2 oo 

= Y,Cn{x;P\q)t^. 



For (7 1 1 they tend to ultraspherical (or Gegenbauer) polynomials: Cn{x;q'^ \ q) — )• C^{x). The 
Gegenbauer polynomials have the generating function 



n=0 

8 Some deeper properties of general orthogonal polynomials 

8.1 True interval of orthogonality 

Consider a system of orthogonal polynomials {pn}- Let pn have zeros Xn,i < Xn,2 < • • • < Xn,n ■ 
Then 

3^1, i ^ 2^1+1,1 > • • • > Xfi^i 4, i^j ^ 00, /o 1 \ 

(o.lj 

Xj^l < < . . . < Xn,n-j+l t < 00. 

The closure I of the interval (^1 ,r]i) is called the true interval of orthogonality of the system 
{pn}- It has the following properties (see [211 p. 11 2]). 

1. I is the smallest closed interval containing all zeros Xn,i- 

2. There is an orthogonality measure for the Pn{x) such that / is the smallest closed interval 
containing the support of /i. 

3. If /i is any orthogonality measure for the Pn{x) and J is a closed interval containing the 
support of fi then I C J. 

8.2 Criteria for bounded support of orthogonality measure 

Recall the three-term recurrence relation (|3.3|) for a system of monic orthogonal polynomi- 
als {pn}- Let be as in (|8.ip . The following theorem gives criteria for the support of an 
orthogonality measure in terms of the behaviour of the coefficients bn, Cn in (I3.3p as n — t- 00. 

Theorem 8.1. 

1. (O p. 109]) If {bn} is bounded and {cn} is unbounded then (^i,7?i) = (— c«,c«). 

2. ([5) Theorem 2.2]) If {bn} and {c„} are bounded then [S,i,r]i] is bounded. 

3. ([5, Theorem 4.5 and p. 121]) If bn ^ b and Cn ^ c {b,c finite) then supp(/i) is bounded 
with at most countably many points outside [6 — 2-^, b + 2y/c], and b^2y/c are limit points 
of snpp{n). 
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Example 8.2. Monic Jacobi polynomials p^'^\x) : 

^ 

" ~ (2n + a + /3)(2n + a + /3 + 2) ^ ' 

_ 4n(n + a)(n + /3)(n + Q + /3) 1 

^ ~ (2n + a + /3 - l)(2n + a + /3)2(2n + a + /3 + 1) ~^ 4 

Hence [6-2^0,6 + 2^] = [-1,1]. 

8.3 Criteria for uniqueness of orthogonality measure 
Put 

p{z) := (Y,\Pn{z)?y {zeC). 



Then < p{z) < oo. Note that p{z) = iff E^olP"(-z)P diverges and that p{z) > iff 
E^o converges. 

Theorem 8.3. ([211 pp. 49-51]) T/ze orthogonality measure is not unique iff p{z) > for all 
z £ C Equivalently, the orthogonality measure is unique iff p{z) = for some z G C. 

In the case of a unique orthogonality measure p, we have p{z) = for z G C\M and p{x) = p{{x}) 
{the mass at x) for x G M, which implies that p{x) iff x is a mass point of p. 

In case of non-uniqueness, for each x G M the largest possible mass of a measure p at x is p{x) 
and there is a measure realizing this mass at x. 

Recall the moments pn := (x",l) = J^x^ dp{x), which are uniquely determined (up to a 
constant factor) by the system {pn}, and also recall the three-term recurrence relation (13. 2p for 
a system of orthonormal polynomials {pn}- 

Theorem 8.4 (Carleman). ( \21\ Theorem 1.10 and pp. 47, 59]) There is a unique orthogonality 
measure for the pn if one of the following two conditions is satisfied. 

oo oo 
n=l n=l 

Example 8.5 (Hermite). 

f°° 2 

P2n = / x'^^e~^ dx = r{n + ^) and log r(n + i) = n log(n + i) + 0(n) as n — oo, 

J — oo 

oo 

so P2n^'^^'^ ~ (n + ^)~2. Hence ^^/"2n''^^"^ ~ °° ' orthogonality measure for the 

n=l 

Hermite polynomials is unique. 
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Example 8.6 (Laguerre). Monic Laguerre polynomials pn satisfy 

xpn{x) = Pn+i{x) + (2n + a + l)pn{x) + n{n + a)pn-i{x). 

oo ^ 

Since > - — tttttt = oo , the orthogonality measure is unique. Also note that 

^ n n + a )V2 

n=0 ^ ^ 

LliQf T{n + a + l) 



hn r(n + l)r(a + l) 

oo 

Since = oo (a > — 1) we conclude once more that the orthogonality measure is unique. 

n=l 

Example 8.7 (Stieltjes-Wigert). Consider the moments /x„ given by the right-hand side of 
Then 

oo oo 
71=1 n=l 

Since the corresponding moment problem is undetermined, the above inequality agrees with 
Theorem I8.4( i). Furthermore, from [12^ (14.27.4)] with q = e~2 we see that the corresponding 
orthonormal polynomials Pnix) = const. Sn{q^x; q) with q = e~2 have = 6^*^(1 — e~2"), by 
which X^^i < oo, in agreement with Theorem 18. 4( ii). 

8.4 Orthogonal polynomials and continued fractions 

Let monic orthogonal polynomials Pn be recursively defined by 

poix) = 1, pi{x) = X-bo, Xpn{x) = Pn+l{x) + bnPn{x) + CnPn-l{x) (fl > 1, Cn> 0). 

Then define monic orthogonal polynomials ^ by 

Po\x) = l, p'l\x)=X-bi, Xp'^\x) = p';^l;^^(x) + bn+lP^nH^) + Cn+lPn-l(.x) (n > 1). 

They are called first associated orthogonal polynomials or numerator polynomials. 
Define 



^IT' ^2(x):= / , Fs{x):= 

x-bo x-bo- x-bo 



x-b2 



and recursively obtain from Fnix) by replacing by H — . This is a 

x-bn 

continued fraction, which can be notated as 

, . 1 |C1 \Cn-2 |c„-l 
Fn[z) = 7 T T T • • • 7 T T • 

z-bo-\z-bi-\ Z- bn-2 - \ z - bn-1 
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Theorem 8.8 (essentially due to Stieltjes). ([5^ Ch. 3, §4]) 



Pn~l{z) (1) ^^.^ _ 1 f Pn{y)-Pn{x) 



^n(2;) = ^^^-r^ and p\'_{y) = — / — d^j.{x). 

Pn{z) fj-o Jm y-x 

Theorem 8.9 (Markov). ([5, Ch. 3, (4.8)]) Suppose that there is a (unique) orthogonality mea- 
sure /i of bounded support for the pn- Let he the true interval of orthogonality. Then 



d^{x) 



z — X 



1 r 

lim Fn{z) = — / 
uniformly on compact subsets o/ C\[^i, ?7i]. 

8.5 Measures in case of non-uniqueness 

Take pn and pn'^ orthonormal: 

Po{x) = l, pi(x) = (x - 6o)/«o, 

XPn{x) = anPn+l{x) + bnPn{x) + an-lPn-l{x) (u > 1), 

p''q\x) = 1, p'i^x) = (x -bi)/ai, 

xp^nHx) = On+ip^j^liix) + bn+iP^Hx) + a„p^j^li{x) (n > 1), 

where > 0. Let = 1, , • • • be the moments for the pn- Suppose that the orthogonahty 
measure for the pn is not unique. Then the possible orthogonality measures are precisely the 
positive measures fi solving the moment problem ()3.4p . The set of these solutions is convex and 
weakly compact. 

We will need the following entire analytic functions. 

oo oo 

.(1) 



A{z) := z^p«(0)p«(z), B{z) := -1 + z ^pW^(0)p„(, 

n=0 n=l 

oo oo 

Ciz) :=l + zY,Pn{0)p^^Uz), D{z)=zY,Pn{Q)Pn{z). 



n=l n=0 

By a Pick function we mean a holomorphic function (j) mapping the open upper half plane into 
the closed upper half plane. Let P denote the set of all Pick functions. In the theorem below 
we will associate with a Pick function (j) a certain measure /U,^. There /it for t G R will mean the 
measure /z,^ with (j) the constant Pick function 2; i— )• t, and /ioo will mean the measure with (j) 
the constant function z 1— 00 (not a Pick function). 

Theorem 8.10 (Nevanlinna, M. Riesz). ( [21^ Theorem 2.12]) Suppose the moment problem 
()3.4p is undetermined. The identity 

df^^jt) _ A(z)0(z) - Cjz) 
t-z B{z)(l){z) - D{z) ^ ' 
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gives a one-to-one correspondence </>—)• /u^ between P U {00} and the set of measures solving the 
moment problem (|3.4p . 

Furthermore the measures Ht it G MU{cxd}) are precisely the extremal elements of the convex 
set, and also precisely the measures fi solving ()3.4p for which the the polynomials are dense in 
L?'{pL). All measures fit are discrete. The mass points of fit are the zeros of the entire function 
tB-D {or of B ift = 00). 

Example 8.11 (Stieltjes-Wigert). The measure which gives in (13. 6p a solution for the moment 
problem associated with special Stieltjes-Wigert polynomials, has a support which is almost 
discrete, but not completely, since is a limit point of the support. Therefore (see the above 
theorem) this measure cannot be extremal. As observed by Christiansen at the end of [6], 
finding explicit extremal measures for this case seems to be completely out of reach. Since the 
measure in (13. 6p is not extremal, the polynomials will not be dense in the corresponding L?' 
space. Christiansen & Koelink [3 Theorem 3.5] give an explicit orthogonal system in this 
space which complements the orthogonal system of Stieltjes-Wigert polynomials to a complete 
orthogonal system. 

9 Orthogonal polynomials in connection with computer algebra 

Undoubtedly, computer algebra is nowadays a powerful tool which many mathematicians and 
physicists use in daily practice for their research, often using wide spectrum computer algebra 
programs like Mathematica or Maple, to which further specialized packages are possibly added. 
This is certainly also the case for research in orthogonal polynomials, in particular when it con- 
cerns special families. Jacobi, Laguerre and Hermite polynomials can be immediately called in 
Mathematica and Maple, while other polynomials in the (g-)Askey scheme can be defined by 
their (g-)hypergeometric series interpretation. Even more geneneral special orthogonal polyno- 
mials can be generated by their three-term recurrence relation. 
Typical kinds of computations being done are: 

1. Checking a symbolic computation on computer which was first done by hand. 

2. Doing a symbolic computation first on computer and then find a hopefully elegant deriva- 
tion which can be written up. 

3. Doing a symbolic computation on computer and then write in the paper something like: 
"By using Mathematica we found . . .". 

4. Checking general theorems, with (hopefully correct) proofs available, for special examples 
by computer algebra. 

5. Formulating general conjectures in interaction with output of symbolic computation for 
special examples. 

6. Trying to find a simple evaluation of a parameter dependent expression by extrapolating 
from outputs for special cases of the parameters. 
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7. Building large collections of formulas, to be made available on the internet, which are fully 
derived by computer algebra, and which can be made adaptive for the user. 

8. Applying full force computer algebra, often using special purpose programs, for obtaining 
massive output which is a priori hopeless to get by hand or to be rewritten into an elegant 
expression. 

While item 8 is common practice in high energy physics, I have little to say about this from 
my own experience. Concerning item 3 there may be a danger that we become lazy, and no 
longer look for an elegant analytic proof when the result was already obtained by computer 
algebra. In particular, many formulas for terminating hypergeometric series can be derived 
much quicker when we recognize them as orthogonal polynomials and use some orthogonality 
argument. 

As an example ot item 1, part of the formulas in the NIST handbook [18j was indeed checked 
by computer algebra. Concerning item 7, it is certainly a challenge for computer algebrists 
how much of a formula database for special functions can be produced purely by computer 
algebra. Current examples are CAOP [1] (maintained by Wolfram Koepf, Kassel) and DDMF 
[8j (maintained by Frederic Chyzak et al. at INRIA). 

The most spectacular success of computer algebra for special functions has been the Zeil- 
berger algorithm, now already more than 20 years old. It is treated in several books: Petkovsek 
et al. [T9j, Koepf [13], Kauers &: Paule [TTj. In particular, p3] contains quite a lot of examples 
of application of this algorithm to special orthogonal polynomials, including the discrete and 
the q-case. 

Various applications of computer algebra to special orthogonal polynomials can be found in 
other chapters of the present volume. 

A very desirable application of computer algebra would be to recognize from a given three- 
term recurrence relation with explict, possibly still parameter dependent coefficients, whether 
it comes from a system of of orthogonal polynomials in the (g-)Askey scheme, and if so, which 
system precisely. A very heuristic algorithm was implemented in the procedure Rec2ortho j20j 
(started by Swarttouw and maintained by the author). It is only up to the level of 2 parameters 
in the Askey scheme. On the other hand Koepf & Schmersau [14] give an algortithm how to 
go back and forth between an explicit eigenvalue equation ()6.ip or (16. 2|) and a corresponding 
three-term recurrence relation with explicit coefficients. 
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